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Abstract 

The neutrino mixing sum rule #12 — #13 cos (5) ~ 01 2 provides a possibility to 
explore the structure of the neutrino mass matrix in the presence of charged 
lepton corrections, since it relates the 1-2 mixing angle from the neutrino mass 
matrix, 1 to observable parameters of the PMNS mixing matrix. The neutrino 
mixing sum rule holds if the charged lepton mixing angles are CKM-like, i.e., small 
and dominated by a 1-2 mixing, and for small 1-3 mixing in the neutrino mass 
matrix. These conditions hold in a wide class of well motivated flavour models. We 
apply this sum rule to present oscillation data, and we investigate the prospects of 
future neutrino facilities for exploring the sum rule by simulating various setups 
for long-baseline reactor and accelerator experiments. As explicit examples, we 
use the sum rule to test the hypotheses of tri-bimaximal and bimaximal neutrino 
mixing, where 0\ 2 i s predicted by sin 2 (^J / 2 ) = 1/3 or 1/2, respectively, although 
the neutrino mixing sum rule can be used to test any prediction for O^i- 
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1 Introduction 



Many attractive classes of models of fermion masses and mixing lead to interesting 
predictions for the neutrino mass matrix m u , such as for instance tri-bimaximal or 
bimaximal mixing. However, the experimentally accessible quantity is the product 
fpMNS = Ve h Vj L . It includes the neutrino mixing matrix V UL , which diagonalizes m u , 
and the charged lepton mixing matrix V £h , which diagonalizes the charged lepton mass 
matrix M e . Often, the essential predictions of flavour models are hidden due to the 
presence of the charged lepton corrections. 

Remarkably, in many cases it can be shown that a combination of the measurable 
PMNS parameters 6*i 2 , 6*i 3 and 5 sums up to the theoretically predicted value for the 1-2 
mixing of the neutrino mass matrix. For example, in an SO (3) family symmetry model 
based on the see-saw mechanism with sequential dominance, predicting tri-bimaximal 
neutrino mixing via vacuum alignment, the neutrino mixing sum rule which is the subject 
of this paper was first observed in the first paper of Ref. pp. In the second and third 
papers of Ref. pQ, it was shown that this neutrino mixing sum rule is not limited to one 
particular model, but applies under very general assumptions, to be specified below. 
For general 1-2 mixing Q\ 2 of the neutrino mass matrix, the neutrino mixing sum rule 
of interest here was given as pQ: 

6 12 -6 13 cos(5) « 6\ 2l (1) 

in the standard PDG parameterization of the PMNS matrix [2] . The specific neutrino 
mixing sum rules for tri-bimaximal [3] and bimaximal [1] neutrino mixing are obtained 
by replacing Q\ 2 by its predictions arcsin(-^) « 35.3° and | = 45°, respectively. 

Let us note at this point that corrections to neutrino mixing angles from the charged 
lepton sector have been addressed in various previous studies. Since some of them are 
sometimes also referred to as "sum rules" , we would like to comment on the differences 
to the neutrino mixing sum rule of Eq. (CQ) [lj. For instance, in many works, it has been 
noticed that charged lepton corrections can make bimaximal neutrino mixing compatible 
with experimental data [SIIEIEIH]- However, in most studies in the literature (e.g. in 
Ref. [S]), CP phases have been ignored. In other works, where complex mixing matrices 
were considered, the connection to the experimentally measurable Dirac CP phase S has 
not been identified [6]. For instance, in Ref. [7j in Eq. (25), the correction is related 
to a phase 0, which is not identical to 5. The introduction of the measurable quantity 
J cp in this equation leads to a sign ambiguity in their "sum rule". We would like to 
remark that various parameterizations of the PMNS matrix are customary, and that 
it is important to specify unequivocally which convention is used. Assuming standard 
PDG parameterization of the PMNS matrix [2] where not stated otherwise, the relations 
of Ref. [8] for bimaximal neutrino mixing are physically inequivalent to the sumrule in 
Eq. ([I]) [lj. Equations for corrections to tri-bimaximal neutrino mixing as well as for 
general 9\ 2 have been considered [HUTU], however the connection to the measurable CP 
phases has not been established. Finally, CKM-like corrections to neutrino mixing angles 
have been considered, and called "sum rules", in Ref. [11J, however CP phases have been 
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ignored. In the following, we will simply refer to the formula of Eq. (TjQ) pQ as the sum 
rule. 

In this paper, after deriving the sum rule [1] , we investigate how well the combination 
of parameters on the left-hand of Eq. (TjQ) can be determined in present and future 
neutrino oscillation facilities, and then compare to the predictions for the right-hand side 
coming from bi-maximal and tri-bimaximal neutrino mixing. Such a study is interesting 
since the sum rule of Eq. (DO) pQ provides a means of exploring the structure of the 
neutrino mass matrix in the presence of charged lepton corrections, and enables whole 
classes of models of neutrino masses and mixings to be tested. However, exploring the 
sum rule requires to measure the currently undetermined mixing angle # 13 as well as the 
CP violating phase 5, which is experimentally challenging, as we shall discuss. 

The outline of the remainder of the paper is as follows. In Sec. |2]we present for the 
first time a parameterization-independent derivation of a family of neutrino mixing sum 
rules, and subsequently show that one of them leads to the sum rule in Eq. (TjQ), using 
the standard PDG parameterization [2] of the PMNS mixing matrix. After this deriva- 
tion of the sum rule, and detailed discussion of the conditions of its validity, we then 
apply it to present oscillation data in Sec. El Sec. H] is devoted to the simulation of fu- 
ture experiments including long-baseline reactor experiments, various second generation 
superbeam setups, a /5-beam, and neutrino factories. We summarize in Sec. El 



2 Derivation of the Sum Rule 

The mixing matrix in the lepton sector, the PMNS matrix t/pMNS; is defined as the 
matrix which appears in the electroweak coupling to the W bosons expressed in terms 
of lepton mass eigenstates. With the mass matrices of charged leptons M e and neutrinos 
m v written as3 

C = -e L M c e R - \v h m v vl + H.c. , (2) 
and performing the transformation from flavour to mass basis by 

K L M e = diag(m e , nfy, m T ), K L rn v Vj L = diag(mi, m 2 , m 3 ), (3) 
the PMNS matrix is given by 

f/pMNS = K L V; t L . (4) 

Here it is assumed implicitly that unphysical phases are removed by field redefinitions, 
and t/pMNS contains one Dirac phase and two Majorana phases. The latter are physical 
only in the case of Majorana neutrinos, for Dirac neutrinos the two Majorana phases 
can be absorbed as well. 

1 Although we have chosen to write a Majorana mass matrix, all relations in the following are 
independent of the Dirac or Majorana nature of neutrino masses. 
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Many attractive classes of models lead to interesting predictions for the neutrino 
mass matrix m v , such as for instance tri-bimaximal [3] or bimaximal [I] mixing where 
VJ h takes the forms 

/ y/2/3 1/VS \ /l/y/2 1/V2 \ 

KU = " VV6 VV3 1/V2 or v£ >u = -1/2 1/2 1/^2 , (5) 
V 1/V6 -1/V3 l/y/2 ) \ 1/2 -1/2 1/^2/ 

respectively, although the sum rule is not necessarily restricted to either of these two 
forms. As mentioned in the introduction such predictions are not directly experimentally 
accessible because of the presence of the charged lepton corrections. However, this 
challenge can be overcome when we make the additional assumption that the charged 
lepton mixing matrix has a CKM-like structure, in the sense that V eh is dominated by 
a 1-2 mixing, i.e. that its elements (V r eL )i3, (K1J23, (K1J31 and (V eL ) 3 2 are very small 
compared to (V eL )ij = 1,2). In the following simplified derivation, we shall take 
these elements to be approximately zero, i.e. 



(KJn (KJ12 

V eL « I (KJ21 (KJ22 I , (6) 
1 



and later on comment on the effect of them being non-zero (see footnote ED- For a 
derivation including these elements, see [TJ] . This situation arises in many generic classes 
of flavour models in the context of unified theories of fundamental interactions, where 
quarks and leptons are joined in representations of the unified gauge symmetries [TUTB"]. 

Under this assumption, it follows directly from Eq. PJ) that (£/pmns)31j (^pmns)32 
and (£7pMNs)33 are independent of V eL , and depend only on the diagonalization matrix 
Vj L of the neutrino mass matrix. This leads to the parameterization-independent sum 
rules which we give in this form for the first time: 

\{Vlhx\ « |(£/pmns)3i| , (7a) 
|«) 32 | « |(C/pmns) 3 2| , (7b) 
\(yl)aa\ « I(^pmns)3 3 |. (7c) 



These innocuous looking relations enable powerful tests of the structure of the neutrino 
mass matrix in the presence of charged lepton corrections. Note that the left-hand sides 
of these relations involve neutrino mixing matrix elements in a particular basis, whereas 
the right-hand sides are basis invariant quantities. This makes sense in the framework 
of a flavour theory which has a preferred basis, the so-called "theory basis" j§ 

2 Also note that models of neutrino masses have a basis-invariant classification. For example, models 
of tri-bimaximal neutrino mixing via Constrained Sequential Dominance (CSD) [14 fall in an invariant 
class of seesaw models, even in the presence of charged lepton corrections, as discussed in [15] . 
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Let us now study the sum rules in the standard PDG parameterization of the PMNS 
matrix (see e.g. [2]), 

C12C13 S12C13 Si 3 e lS 

UpMNS = I -C23S12 - Si 3 S 2 3Ci 2 e l5 C23C12 - Si 3 S 2 3S 12 e l8 S23C13 I P M aj , (8) 

S23S12 - si 3 c 23 c 12 e lS -S23C12 - si 3 c 23 si2e l5 c 23 ci 3 

which is used in most analyses of neutrino oscillation experiments. Here S is the so- 
called Dirac CP violating phase which is in principle measurable in neutrino oscillation 

.01.02 

experiments, and P^aj = diag(e : 2 , e 1 2 , 0) contains the Majorana phases ai, a 2 . In the 
following we will use this standard parameterization (including additional phases) also 
for Vj L and denote the corresponding mixing angles by 9^, while the mixing angles 9ij 
without superscript refer to the PMNS matrix. 

In addition to the assumption that V eL is of the form of Eq. we will now assume 
that the 1-3 mixing in the neutrino mass matrix is negligible, 

MDi8«0. (9) 

Many textures for the neutrino mass matrix fulfill this relation exactly, for example the 
cases of bimaximal and tri-bimaximal mixing, although the assumption in Eq. @ is 
more general. Using the assumption © in the sum rule of Eq. (I7a|) one obtains 

S 23 S 12 ~ |s 2 3Sl2 - Si 3 C23Ci 2 e 1,5 | « S 23 S 12 ~ S 13 C 23 C 12 COs(<5) , (10) 

where the last step holds to leading order in si 3 . Furthermore, Eq. (ITcl) together with 
Eq. (jUJ) implies 

e v 23 « 23 + o(e 2 l3 ) . (ii) 

Using this relation in Eq. (flOl) leads to the sum rule 

12 -0 13 cot(0 23 )cos(5)«^ 2 , (12) 



which holds up to first order in # 13 . Hence, we have obtained an approximate expression 
for the (in general unobservable) mixing angle d\ 2 in terms of directly measurable pa- 
rameters of the PMNS matrix. This sum rule can be used to test a bimaximal [Q\ 2 = ?) 
or tri-bimaximal (6^ 2 = arcsin(^)) structure of the neutrino mass matrix, but may as 
well be applied for a different pattern beyond these two examplesjf] In the following we 
will specialise our discussion to models predicting maximal 2-3 mixing in the neutrino 
mass matrix, Q 23 = \- This includes of course the cases of bimaximal and tri-bimaximal 
mixing. With Eq. (iTTj) this leads to the sum rule of Eq. ([1]) [1] . 

3 We would like to remark at this point that the sum rule holds at low energy, where the neutrino 
oscillation experiments are performed. Therefore, if theory predictions arise at high energies like the 
GUT scale, their renormalization group evolution has to be taken into account. In seesaw models, the 
running can be calculated conveniently using the software package REAP [IT)] . 
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As a side remark we mention that under the above conditions also a simple rela- 
tion for |(Z7pmns)i3| can be obtained, |(£/pmns)i3| ~ | (V^ L )i 2 1 KVjJ 23 !- m tne stan dard 
parameterization it yields 

9 13 ^s 23 9 e 12 ^j= 2 9 e 12 , (13) 

where | (V^ L ) 12 1 = sm(9l 2 ) m 9f 2 , and in the last step we have approximated 9 23 ~ f. 
Hence, #13 is related to the 1-2 mixing in the charged lepton mass matrix]^ This relation 
has been noticed by many authors, e.g. Refs. [IJEllQdIBllHlIISIIIlIIB] , an d it can provide 
additional hints on the underlying theory of flavour (see for example Ref. pQ). Here 
we will not explore this relation further but focus on the sum rule ([T]). Note, however, 
that under the assumptions and the only mixing parameters of the model are 
dl 2 ,9i 2 ,9 23 , and through the relations (TT3"|) . (TT2"]) . (TTT]) all of them can be expressed in 
terms of measurable PMNS parameters. 

Finally, we mention that under the above assumptions Eq. ( jTTj) can be used to test 
predictions for 9 23 . This is complementary to the application of the sum rule for 9" 2 , 
and a precise determination of 9 23 will allow for an additional test of predictions for 
the neutrino mass matrix^ For the examples of tri-bimaximal and bimaximal neutrino 
mixing, one can test experimentally the prediction #23 ~ f , in addition to the verification 
of the corresponding sum rule for 6*12. Prospects for the measurement of deviations from 
maximal 2-3 mixing have been discussed e.g. in Ref. [PT] . 



3 The Sum Rule and Present Oscillation Data 

In this section we show that already with present global data from neutrino oscilla- 
tion experiments the sum rule can be used to test the hypotheses of bimaximal or 
tri-bimaximal mixing in the neutrino mass matrix. Using the results from the global 
analysis of Ref. [18] a fit is performed under the assumption of the sum rule Eq. (CD) 
with the constraints 9\ 2 = 45° for bimaximal mixing or 9 V X2 = arcsin(^) « 35.3° for 
tri-bimaximal mixing. (See also Ref. [19] for similar considerations.) 

Present data implies that 9 12 is significantly smaller than 45°, with the upper limit of 
39.2° at 3cr dominated by the SNO solar neutrino experiment [20]. Hence, to reconcile 
the value 9\ 2 = 45° for bimaximal mixing one needs a relatively large value of (9 13 and 
cos(<5) ~ —1. These expectations are confirmed by the fit, as visible in the left panel 
of Fig. [U We obtain Ax^in ~ 6.14 with respect to the unconstrained best fit point, 
and hence, allowed regions appear only at 99% and 99.73% CL. Fitting the sum rule for 

4 Note that in the derivation of the sum rule Eq. (TTJ)) it was not necessary to assume that 9\ 2 is 
small. This requirement follows only a posteriori from Eq. (|13p and the fact that #13 has to be small 
from data. 

If the assumption of Eq. ([6]) is relaxed and one allows for a small (but non-zero) 2-3 mixing in 
the charged lepton mixing matrix, 6*| 3 -C 1, there will be a correction of order 6*2 3 to Eq. (TTT|) . which 
has to be taken into account when drawing conclusions on #23 from a measurement of 623. It can be 
shown p, however, that this correction does not affect the #12 sum rule (fT2"|) to leading order. 
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Figure 1 : Fit to the present global data on neutrino oscillations under the assumption of the sum rule 
and bimaximal (left) and tri-bimaximal (right) mixing in the neutrino mass matrix. Allowed regions 
are shown at 90%, 95%, 99%, 99.73% CL (2 dof) with respect to the unconstrained best fit point. The 
dashed lines correspond to the upper bound on #13 at 3ct. 

bimaximal mixing requires that both, # 12 and #13, are pushed towards their upper limits 
which leads to the increase of y 2 mentioned above. We conclude that already present 
data disfavours bimaximal neutrino mixing under the assumption of the sum rule at 
more than 2a. Hence, either the hypothesis of 9" 2 = 45° has to be discarded, or some of 
the approximations needed for the sum rule are not justified, for instance the charged 
lepton corrections cannot be of CKM type as assumed in Eq. ©. On the other hand, 
if the fit is accepted, the sum rule for bimaximal mixing predicts that #13 is close to its 
present bound and S ~ ir. 

The right panel of Fig.fUshows the result for the sum rule with tri-bimaximal neutrino 
mixing. In this case the fit is fully consistent with the data and a best fit point at the 
same \ 2 as the unconstrained fit is found. This follows since the best fit point # 12 = 33.2° 
is close to the tri-bimaximal mixing value. Indeed, for small values of 8 13 , say less than 
2°, the sum rule is satisfied within current experimental errors, for all values cos (5). On 
the other hand, the sum rule leads to a strengthening of the bound on #13 in the regions 
where cos(<5) 7^ 0. 



4 The Sum Rule and Sensitivities of Future Exper- 
iments 

In this section we explore the ability of future experiments to constrain the parameter 
combination Qyi — 6>i 3 cos(5), appearing on the left-hand side of the sum rule in Eq. (Tj[|), 
in order to obtain information on 6^ 2 , and so enable a comparison with the predicted 
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accuracy on 9u at 3 a 


pieseill Lid Id 


A Q° 


SK-Gd, 5 years 


4.0° 


LENA © Frejus, 44 kt, 5 y 


2.0° 


SPMIN @ 60 km, 200 GW kt y 


1.1° 



Table 1: Accuracy on 812 from present data and three options for future LBL reactor experiments. 

values of 9^ 2 coming from particular flavour models. Obviously, in order to do this, the 
errors on 812 as well as on #13 cos(<5) should be as small as possible. 

Let us first discuss the prospects to improve the accuracy on 9%2, which is 4.9° 
at 3ct from present data [IE], dominated by the SNO solar neutrino experiment [20J. 
Significant improvement on 9 12 can be obtained by long-baseline (LBL) reactor neutrino 
experiments, similar to the KamLAND experiment [21J. A realistic possibility is that 
the Super-K experiment will be doped with Gadolinium (SK-Gd) [22] • This will allow 
for a very efficient detection of reactor u e , and by observing the neutrino flux from 
the surrounding reactors in Japan a precise determination of the "solar" oscillation 
parameters will be obtained [23]. Following the analysis of Ref. [21], after 5 years of 
data taking an accuracy of 4.0° can be obtained for # 12 at 3 ex. Another interesting 
option would be a big scintillator detector such as LENA. In Ref. [21] the possibility of 
a 44 kt detector installed in the Frejus underground laboratory has been considered. By 
the observation of the reactor neutrino flux from nearby reactors in France an accuracy 
of 2.0° can be obtained for 9\2 at 3 o after 5 years of data taking. 

Probably the best way to measure #12 would be a dedicated reactor experiment 
with only one reactor site at a distance around 60 km [2H1I2E], where the first survival 
probability minimum would be right in the middle of the reactor event rate spectrum. 
This has been named "SPMIN experiment" in Ref. [26]. The obtainable accuracy in this 
type of experiment, as in all reactor experiments, is a balance between statistical and 
systematical errors. The former call for large detectors and powerful reactors, whereas 
the latter require great experimental skill and a careful design. For illustration we 
consider here a rather "big'' setup corresponding to an exposure of a liquid scintillator 
detector of 200 GW kt yO The estimated accuracy at 3 a of such an experiment to 
9\2 is 0.7° from statistical errors only, and 1.1° if various systematical effects are taken 
into account. These numbers have been obtained by applying a similar analysis for 
the SPMIN experiment as in Ref. [21], where also a detailed description of the various 
systematics can be found. At such large exposures systematics have a big impact on the 
accuracy, but it seems difficult to improve the systematics in a very large kiloton sized 
detector. 

The reason why neither SK-Gd nor LENA at Frejus can compete with a dedicated 
SPMIN experiment is that many nuclear reactors at various distances contribute which 
washes out the oscillation signature to some extent. Let us note that also future solar 

6 For comparison, typical nuclear power plants have a thermal power output of order 10 GW, and 
the KamLAND experiment has a total mass of about 1 kt. 
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experiments, even with a 1% measurement of the pp-flux, cannot compete with an 
SPMIN reactor experiment [26]. The prospects of the 812 measurement are summarized 
in Tab. HJ In the following combined analysis with LBL accelerator experiments we will 
consider as reference values the 3 a accuracies of 4.0°, 2.0° and 1.1° obtainable at SK-Gd, 
LENA at Frejus, and a 200GWkty SPMIN experiment, respectively. 

Next we turn to the sensitivity of LBL accelerator experiments to the combination 
of physical parameters appearing on the left-hand side of the sum rule in Eq. (TjQ) . These 
experiments are sensitive to 5 and #13 but have nearly no sensitivity to 812- Therefore, 
we will use the input on 812 from LBL reactor experiments as described above and 
perform a combined reactor plus accelerator analysis. We follow the general analysis 
procedure as described in Ref. with the difference that we now project onto the 
direction 8yi — #13 cos(<5) in the parameter space. Thus, the obtained results do include 
the errors and the correlations on #13 and 8 as well as the errors on 6*12, Am^, 6*23; Am 2 ^ 
and the matter density. Especially the correlation between #13 and 8 is crucial, since 
the relevant oscillation probability contains terms which go as 

8 13 sm(5) and 6>i 3 cos(5). (14) 

However, the ///independence of these two terms is different and hence experiments 
covering different L/E'-ranges may have very different sensitivities to #i 3 cos(<5). For 
these reasons the accuracy on the combination 13 cos(<5) may be very different from the 
accuracy individually obtained on 813 and 8. Therefore, a proper treatment and inclusion 
of the correlation between 9 is and 8 is mandatory to obtain meaningful results. 

For the experiments discussed in the following the sensitivity to either 813 sm(5) 
or #i3 cos(5) is dominated by the data from the appearance channels. There are two 
main reasons for this: The sin(<5) term can only appear in off-diagonal transitions, i.e. 
appearance channels, because it is manifestly CP violating. Secondly, a possible cos(<5) 
contribution in the disappearance channels is always suppressed with respect to the 
leading #23 effect and hence plays no statistically significant role. Only for the very 
largest values of #13 there is contribution of the disappearance channels, but it is still 
very small. We will not discuss the possible impact of short baseline reactor experiments 
which are designed to determine 6 I 13 . The reason is, that all the experiments discussed 
in the following have a superior sensitivity to # 13 on their own. 

The calculations are performed with the GLoBES software package [2E]- For the 
input values of the oscillation parameters we use Am^ = 2.5 x 10~ 3 eV 2 , Am^ = 
8 x 1CT 8 eV 2 , #23 = 45°, and 812 = 33.2°. We consider six examples for future experi- 
ments. Their main characteristics are summarized in Tab. [21 They include three second 
generation superbeam experiments, SPL - a CERN based experiment with a Mt size 
water Cerenkov detector at Frejus [2H], T2HK, the second stage of the Japanese T2K 
project [29J (see also Ref. [27j). and WBB - a wide-band beam with a very long baseline 
as discussed in the US [30,31j. Furthermore, we consider an advanced /3-beam setup 
BB350 as described in Ref. [32], with a relativistic 7-factor of the decaying 18 Ne and 
6 He ions of 350. All these experiment are planed to employ a large water Cerenkov 
detector with a fiducial mass in the range 300 — 440 kt. Note, that the setup labeled 
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Setup 


Ref. 


Baseline 


Detector 


Beam 


SPL 


m 


130 km 


440 kt WC 


4 MW superbeam, 2 y (v) + 8 y (p) 


T2HK 


[29\ 


295 km 


440 kt WC 


4 MW superbeam, 2 y (v) + 8 y (p) 


WBB 


w\ 


1300 km 


300 kt WC 


1.5 MW superbeam, 5 y [v) + 5 y (v) 


BB350 


[32 


730 km 


440 kt WC 


5x1.1- 10 18 18 Ne + 5 x 2.9 • 10 18 6 He 


NFC 


m 


4000 km 


50 kt MID 


50 GeV, 4 x 10 21 yT + 4 x 10 21 


NFO 


[33] 


4000+7500 km 


2x50 kt MID* 


20 GeV, 4 x 10 21 fT + 4 x 10 21 /x+ 



Table 2: Summary of the six future LBL accelerator experiments considered in this study. WC 
stands for water Cerenkov detector and all masses for this technology are fiducial masses. MID denotes 
a magnetized iron calorimeter, whereas MID* denotes an improved version thereof. In the column 
"Beam" we give for BB350 the total number of useful ion decays, and for NFC, NFO the energy of the 
stored muons and the total number of useful muon decays. For more details see the text. 



WBB assumes an operational time per solar year of 1.7 • 10 7 s instead of the usual 10 7 
s. The two neutrino factory setups considered here, NFC and NFO are taken from [33J . 
NFC is what we call conservative, in the sense that it employs only one magnetized iron 
detector (MID) with the canonical properties regarding muon detection threshold and 
background rejection [SUEZ]. NFO is an optimized version, which uses two identical 
detectors at two baselines of 4000 km and 7500 km, the latter being the so-called magic 
baseline [35]. The second difference is that the detector is now an improved MID*, 
which has a lower muon detection threshold but somewhat larger backgrounds, for de- 
tails see [33]. The lower threshold allows to reduce the muon energy from 50 GeV to 20 
GeV. 

Figs. [5] and [3] show the results for the 3cr allowed interval in 8 12 — 9i 3 cos(5) as 
a function of the true value of 5 assuming sin 2 2^13 = 10 _1 and 10~ 2 , respectively, 
from the considered experimental setups. This allowed interval can be compared with 
theoretical predictions for Q\ 2 . We illustrate in the figures the cases of bimaximal and 
tri-bimaximal mixing by the horizontal lines, but of course any prediction for Q v l2 can be 
confronted with the outcome of the experiments. Since we have used as true value for 
Qyi the present best fit point of 33.2°, bimaximal mixing {Q\ 2 = 45°) can be obtained 
only for large values of #13 and 5 ~ 180°, in agreement with the discussion in Sec. [3l For 
larger (smaller) true values of 9\ 2 , the bands and islands in Figs. [2] and [3] are shifted up 
(down) correspondingly. 

All experiments shown in Figs. [2] and [3] have good sensitivity to 612 — #i3Cos(5). 
In many cases only some specific values of the CP phase 5 are consistent with a given 
prediction for 9" 2) which illustrates the power of the sum rule. An interesting observation 
is that the presence of the mass hierarchy degenerate solutions (dashed lines) limits the 
usefulness of SPL and T2HK severely. In these experiments the matter effect is small 
because of the relatively short baseline. This implies that the mass hierarchy degenerate 
solution cannot be resolved. Furthermore, the degenerate solution appears at a similar 
value of 6> 13 but at a fake CP phase close to n — 5 [36J. This changes the sign of the 
term # 13 cos(5), which explains the shape of the dashed curves in the figures. Because of 
this degeneracy an ambiguity appears when the sum rule is applied for SPL and T2HK, 
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Fi gure 2: The 3<7 allowed interval for the parameter combination 6\2 — #13 cos(<5) appearing on the 
left-hand side of the sum rule in Eq. (TJJ), as a function of the true value of S for sin 2 26*13 = 10 _1 from 
various LBL experiments. The dashed lines correspond to the sgn(Am§ 1 ) degenerate solution. The 
colors indicate different errors on 612'- blue - 4.0°, red - 2.0° and green - 1.1° (at 3a). A true value 
sin 2 #12 = 0.3 (6*12 = 33.2°) has been assumed. The horizontal lines show the sum rule predictions 
corresponding to bimaximal and tri-bimaximal neutrino mixing. 



which significantly limits the possibility to distinguish between various predictions for 
6" 2 , especially for large values of #13 as visible in Fig. [21 A solution to this problem 
could be the information provided by atmospheric neutrinos in the Mt size detectors 
used in these experiments [37J (which is not included here). For the other experiments 
the problem of the degeneracy is absent, since the mass hierarchy degeneracy can be 
resolved (at sufficiently large #13) thanks to the longer baselines. 

The performance of all experiments is summarized also in Fig. HJ which shows the 
obtainable 3 a accuracy for the combination of parameters # 12 — # 13 cos(<5), appearing 
on the left-hand side of the sum rule in Eq. (|T|), as a function of the true value of 5 for 
the two cases sin 2 2# 13 = 1CT 2 and 1CT 1 . This figure shows that it will be possible to 
discriminate between models whose predictions for Q\ 2 differ by a few degrees. For the 
large value of #13 assumed in Figs. [2] and H] (right), sin 2 2# 13 = 1CT 1 , the total uncertainty 
is dominated by the term # 13 cos(5) in the sum rule, and a modest improvement of the 
current error on # 12 will be enough for exploring the sum rule. The accuracy depends 
significantly on the true value of 5. Obviously the impact of the term # 13 cos(<5) is 
larger for |cos(5)| ~ 1. For smaller values of 813 the accuracy on 612 becomes more 
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Figure 3: Same as Fig. fjbut for sin 2 26 13 = 1(T 2 . 



important, the overall sensitivity is dominated by the LBL reactor measurement, and 
the dependence on 5 is weaker \j 

It follows from Figs. El El and H] that NFO has the best performance for all values 
of 6* 13 , making this the machine of choice for testing the sum rule. NFC compares well 
to BB350 for this measurement, whereas the performance on 5 and 6* 13 individually 
is much worse for NFC than for BB350. The reason for this behaviour is that an 
experiment whose events are centered around the first oscillation maximum like a /3-beam 
or superbeam is mainly sensitive to the # 13 sin(r5) term. A neutrino factory, however, 
gets most of its events above the first oscillation maximum and thus is much more 
sensitive to the 6* 13 cos(5) term. This explains also the relatively good performance of 
the WBB visible from the right panel of Fig. HI where WBB performs second only to 
NFO. For such large values of #13 (sin 2 29^ = 10 _1 ) spectral information far beyond the 
first oscillation maximum can be explored efficiently, which is important for constraining 
6>i3 cos(<5). 

For the somewhat smaller value of # 13 , sin 2 2# 13 = 10~ 2 , the performances of NFO, 
NFC, and BB350 become rather similar, whereas the accuracies obtainable at super- 
beams depend still to some extent on the true value of 5, see Fig. 0] (left). Note that in 

7 As visible in Fig. [3l for a few values of 5 for T2HK the allowed region of 6*12 — 6*13 cos(<5) consists 
of two disconnected intervals even for fixed neutrino mass hierarchy, because of the so-called intrinsic 
degeneracy. This explains the "turn over" of the T2HK lines at some values of 8 in Fig. H] (left). 
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Figure 4: The 3 a error in degrees for the combination of parameters #12 — #13 cos(<5), appearing on the 
left-hand side of the sum rule in Eq. (fTJ), as a function of the true value of S for sin 2 2#i3 = 1CP 2 (left) 
and sin 2 2^i3 = 1CP 1 (right). The different colored lines are for different experiments as given in the 
legend. The sgn(Am 2 1 ) degenerate solution has been omitted. Although sin 2 612 = 0.3 (6*12 = 33.2°) 
has been used as true value, the results are practically independent of this assumption. The error on 
6»i2 is 1.1° at 3 a. 



this figure the most optimistic accuracy on 612 from an SPMIN reactor experiment has 
been assumed, and that the mass hierarchy degeneracy has not been taken into account. 
Indeed, decreasing the true value of #13 for all setups besides the neutrino factory one, at 
some point the mass hierarchy degenerate solution kicks in and introduces an ambiguity 
in the allowed interval for 6*12 — #13 cos(5), compare also Fig. [3l 

For a given model prediction of the neutrino mixing angle 6" 2 , the sum rule in Eq. 
may be expressed as a prediction for the physical solar mixing angle as a function of 
the CP violating Dirac oscillation phase 5. Fig. [5] shows the sum rule prediction for 
the PMNS parameter # 12 corresponding to tri-bimaximal mixing in the neutrino mass 
matrix, Q\ 2 = arcsin(^) rs 35.3°, i.e., 

U » 35.3° + i3 cos(<J) . (15) 

In the figure we have simulated data for the NFO setup for different true values of #13 
and 5 and used Eq. ( TT51) to calculate the resulting 3 a range for the predicted 9i 2 . This 
result can be compared with the outcome of a separate measurement of #12 (for example 
in a reactor experiment) to test whether the hypothesis of tri-bimaximal neutrino mixing 
is compatible with the assumptions leading to the sum rule. 
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Figure 5: The prediction for 612 from the sum rule with tri-bimaximal neutrino mixing, as given in 
Eq. (|15p . The different curves show the predicted 3 a interval for #12 following from a NFO measurement 
of #13 cos(5) as a function of the true values of 8 and #13. 

5 Summary and Conclusions 

In this work we have considered the sum rule in Eq. (JTJ, and in particular how well 
the combination of parameters 6 , 12 — #13 cos(5), which appears on the left-hand side, can 
be measured in oscillation experiments. This is important, since the sum rule follows 
from quite general assumptions which are satisfied in a wide class of flavour models. 
Moreover, particular such flavour models make definite predictions for #i 2 , an d the sum 
rule then enables these models to be tested. 

We have derived the sum rule, starting from a parameterization independent set 
of sum rules, which follow from certain well defined assumptions about the nature of 
charged lepton and neutrino mixings. We then expressed the sum rule in terms of 
the standard PMNS mixing parameters (see e.g. [2]) commonly used in presenting the 
results of neutrino oscillation experiments. One way to view the sum rule is to consider 
the charged lepton corrections to the neutrino mixing angle #i 2 predicted from theory, 
leading to the physical solar neutrino mixing angle #12. Then, under certain assumptions, 
the charged lepton correction turns out to only depend on the physical combination 
#13 cos(5). To be precise, the sum rule in Eq. ([T]) holds up to first order in #13 under the 
following assumptions: 

(a) The charged lepton mixing matrix is CKM-like, i.e., dominated by the 1-2 mixing 
angle, see Eq. 

(b) The 1-3 element of the neutrino mixing matrix is negligible, #i 3 ~ 0. 
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(c) The 2-3 mixing in the neutrino mass matrix is maximal, ~ 45°, which under 
the previous two assumptions is equivalent to #23 ~ 45° + 0(9f 3 ). 

If condition (c) is not satisfied and 823 turns out to be non- maximal the generalized 
sum rule Eq. ( 1121) has to be used. This would not change the reasoning of this paper, 
since 8\ 2 still can be expressed only in terms of measurable quantities, involving now 
also 623 ■ Hence, condition (c) has been adopted here just for simplicity. Condition (a) 
holds for a large class of models. For example, in many GUT models the charged lepton 
mixing matrix is related to the quark mixing matrix since quarks and leptons are joined 
in representations of the unified gauge symmetries. On the other hand, conditions (b) 
and (c) are rather typical for flavour models in the neutrino sector. In particular, the 
popular examples of bimaximal and tri-bimaximal mixing fulfill conditions (b) and (c) 
exactly. 

We have demonstrated the usefulness of the sum rule by imposing it as a constraint in 
a fit to present global data from neutrino oscillation experiments under the assumptions 
of bimaximal and tri-bimaximal neutrino mixing. This analysis shows that under the 
condition (a) bimaximal neutrino mixing is disfavoured at about 2cr by present data 
with respect to tri-bimaximal mixing, which is perfectly compatible with the data. If 
the fit for bimaximal mixing is accepted the sum rule predicts that 613 is close to its 
present bound and 5 ~ tt. 

In the main part of the paper we have concentrated on the sum rule in the context 
of future high precision neutrino oscillation experiments. We have considered long- 
baseline reactor experiments for a precise measurement of 812 , as well as six examples for 
advanced long-baseline accelerator experiments to constrain the parameter combination 
6*i3 cos(<5) appearing in the sum rule. These setups include three options for second 
generation superbeam experiments, a /3-beam, and two examples for a neutrino factory. 
It is shown that most of these experiments will allow for a rather precise testing of the 
sum rule, and 8 V 12 can be inferred within an accuracy of few degrees, where the precise 
value shows some dependence on the true values of #13 and 5. For sin 2 2^13 < 10~ 2 the 
accuracy is dominated by the error on 812, whereas for large values of 813 the precision on 
the term #13 cos(5) dominates. Obviously its impact is larger for | cos(5)| ~ 1. Because of 
the appearance of cos(5) experiments operating not only at the first oscillation maximum 
(where there is good sensitivity to sin(5)) are well suited for this kind of measurement, for 
example a neutrino factory or a wide-band superbeam. Another interesting observation 
is that the mass hierarchy degeneracy plays an important role. Since this degeneracy 
introduces an ambiguity in the CP phase 5 its appearance significantly reduces the 
information on 8\ 2 which can be extracted from the sum rule. 

To conclude, the neutrino mixing sum rule considered in this work is a convenient 
tool to explore the structure of the neutrino mass matrix in the presence of charged 
lepton corrections, and to test whole classes of models of neutrino masses and mixings. 
Already applied to present data it is possible to obtain non-trivial statements, whereas 
with future high precision oscillation experiments a rather accurate testing of models 
will become possible in the framework of the sum rule. 
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